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^ We study logarithmically enhanced electromagnetic corrections to the decay rate in the high 
dilepton invariant mass region as well as corrections to the forward backward asymmetry (FBA) 
of the inclusive rare decay B Xsi~^i~ . As expected, the relative effect of these corrections 
^ [ in the high dilepton mass region is around —8% for the muonic final state and therefore much 
' larger than in the low dilepton mass region. 
We also present a complete phenomenological analysis, to improved NNLO accuracy, of the 
dilepton mass spectrum and the FBA integrated in the low dilepton mass region, including a 
new approach to the zero of the FBA. The latter represents one of the most precise predictions in 
flavour physics with a theoretical uncertainty of order 5%. We find (qq)^^ = (3.50±0.12)GeV^. 
For the high dilepton invariant mass region, we have B{B Xs/i/i)high = (2.40lo;g2) x 10^^. 
The dominant uncertainty is due to the 1/mb corrections and can be significantly reduced in 
the future. For the low dilepton invariant mass region, we confirm previous results up to small 
corrections. 
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1 Introduction 



At the beginning of the LHC era, the search for physics beyond the Standard Model (SM) 
is the main focus of particle physics. In principle, there arc two ways to search for possible 
new degrees of freedom. At the high-energy frontier one tries to produce those new degrees of 
freedom directly, while at the high-precision frontier one analyses the indirect virtual effects of 
such new particles. It is a matter of fact that high-precision measurements allow to analyse 
new physics scales presently not accessible at direct collider experiments. 

There are theoretical arguments like the hierarchy problem which let us expect new physics 
at the electroweak scale. However, the indirect constraints on new physics by the present 
flavour data indicate a much higher new physics scale when such new effects are naturally 
parametrized by higher-dimensional operators. Thus, if there is new physics at the electroweak 
scale, then its flavour structure has to be highly non-trivial and the experimental measurement 
of flavour- violating couplings is mandatory. This 'flavour problem' has to be solved by any new 
physics scenario at the electroweak scale. Moreover, the present electroweak data also indicate 
a slightly higher new physics scale ('little hierarchy problem'). 

Among inclusive flavour-changing neutral current (FCNC) processes (for a review see [1,2]), 
the inclusive B X^i^i^ decay presents an important test of the SM, complementary to the 
inclusive B Xs^ decay. It is particularly attractive because of kinematic observables such 
as the dilepton invariant mass spectrum and the forward-backward asymmetry (FBA). These 
observables are dominated by perturbative contributions if the cc resonances that show up as 
large peaks in the dilepton invariant mass spectrum are removed by appropriate kinematic cuts. 
In the so-called 'perturbative g^-windows', namely in the low-dilepton-mass region 1 GeV^ < 

= mjf^ < 6GeV^, and also in the high-dilepton-mass region with > 14.4 GeV^, theoretical 
predictions for the invariant mass spectrum are dominated by the perturbative contributions, 
and a theoretical precision of order 10% is in principle possible. 

The integrated branching ratio has been measured by both Belle [3] and BaBar [4] based on 
a sample of 152 x 10^ and 89 x 10^ BB events respectively. In the low-dilepton invariant mass 
region, 1 GeV^ < g^ < 6 GeV^, the experimental results read 

r (1.493 ±0.504l°:il) X 10-6 (Belle) 
B{B X,£+£-)iow = I (1.8 ± 0.7 ± 0.5) x 10-^ (BaBar) (1) 

I (1.60 ± 0.50) X 10-6 (Average) 

Measurements for the high-g^ region, 14.4 GeV^ < g^ < 25 GeV^, are also available: 
B{B ^ X,£+r )high 

By the end of the present B factories an experimental accuracy of 15% is finally expected. 

The recently calculated NNLL QCD contributions [5-14] have significantly improved the 
sensitivity of the inclusive B —>■ Xst^t~ decay in testing extensions of the SM in the sector 



(0.4I8 ± CllTi^^y X 10-6 (Belle) 

(0.5 ± 0.251H?) X 10-6 (BaBar) (2) 

(0.44 ± 0.12) X 10-6 (Average) 
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of flavour dynamics, in particular, the value of the dilepton invariant mass for which the 
differential FBA vanishes is one of the most precise predictions in flavour physics with a theo- 
retical uncertainty of order 5%. This well corresponds to the expected experimental sensitivity 
of 4 — 6% at the proposed Super- i? factories [15-17]. 

Also non-perturbative corrections scaling with l/ml, l/ml, or 1/m^ [18-24] have to be 
taken into account. Moreover, factorizable long-distance contributions away from the resonance 
peaks are important; here using the Kriiger-Sehgal approach [25] avoids the problem of double- 
counting. 

In the high-g^ region, one encounters the breakdown of the heavy-mass expansion at the 
endpoint; while the partonic contribution vanishes section in the end-point, the l/ml ^^'^ l/ml 
corrections tend towards a non-zero value. In contrast to the endpoint region of the photon 
energy spectrum in the B — > Xg^ decay, no partial all-order resummation into a shape function 
is possible here. However, for an integrated high-g'^ spectrum an effective expansion is found 
in inverse powers of mf^ = m?, x (1 — y/smin) rather than mi, [26,27]. The expansion converges 
less rapidly, depending on the lower dilepton mass cut Smin [11]. Recently it was suggested [24] 
that the large theoretical uncertainties can be reduced by normalizing the B — > Xs^'^^~ decay 
rate to the semileptonic B X„£z/ decay rate with the same cut. 

A hadronic invariant- mass cut is imposed in the present experiments (Babar: mx < 1.8 GeV, 
Belle: mx < 2.0 GcV) in order to eliminate the background such as h ^ c (— se^v)e~i' ~ 
b se^e~ + missing energy. The high-dilepton mass region is not affected by this cut, but 
in the low-dilepton mass region the kinematics with a jet-like Xs and m\ < msAqco implies 
the relevance of the shape function. A recent SCET analysis shows that using the universality 
of jet and shape functions the 10 — 30% reduction of the dilepton mass spectrum can be 
accurately computed using the B — > shape function. Nevertheless effects of subleading 
shape functions lead to an additional uncertainty of 5%. It is also shown that the zero of the 
FBA is not affected by the cut [28,29]. We emphasize that, as in all previous phenomenological 
analyses of the inclusive B — > mode, the predictions in this paper are given under the 

assumption that there is no cut in the hadronic mass region. 

Finally, as in the B — > case, there are unknown subleading nonperturbative corrections 
of order 0{asA./mf,) which may be estimated by an additional uncertainty of order 5%. 

Recently, further refinements were presented such as the NLO QED two-loop corrections 
to the Wilson coefficients whose size is of order 2% [12]. Furthermore, it was shown that 
in the QED one-loop corrections to matrix elements large coUinear logarithms of the form 
log(m^/mfgp^.Qj^) survive integration if only a restricted part of the dilepton mass spectrum is 
considered. This adds another contribution of order +2% in the low-g^ region for B{B 
Xs^'^^") [30]. For B{B XgC^e"), in the current BaBar and Belle setups, the logarithm of 
the lepton mass gets replaced by angular-cut parameters and the integrated branching ratio for 
the electrons is expected to be close to that for the muons. 

In the present manuscript, we calculate these corrections due to large coUinear logarithms 
also for the dilepton mass spectrum in the high-g^ window and for the FBA for the first time. 
In the high-g^ region these corrections are much larger than in the low-g^ one; we already 
anticipate that in the muon channel the former account for —8% of the rate while the latter 
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are only about +2%. We also update the theoretical prediction of the observables using several 
improvements offered recently in the literature among which the two-loop matrix element of the 
semileptonic operator Pg and the 1/ml corrections are the most important to date. We present 
a new phenomenological analysis of the dilepton mass spectrum in the low- and high-g^ region 
and of the FBA. In particular, a new approach to the zero of the FBA is proposed. Moreover, 
we explore a strategy to reduce the uncertainties due to the l/mi, corrections by normalizing 
with the semileptonic B X^iu decay rate with the same cut, as was recently proposed 
in the literature [24]. We already anticipate that with this normalisation the dilepton mass 
spectrum, integrated over the high-g^ region, receives only a 3% shift in the central value by 
the nonperturbative corrections. 

Finally, we derive model-independent formulae for non-SM values of the high-scale Wilson 
coefficients. 

The present manuscript is organised as follows. In Section [5] we present the master formulae 
for the calculation of the observables; in particular we discuss the normalization of the branching 
ratio and forward backward asymmetry and explain the details of the perturbative expansion. 
In Section [3] we present the novel analytical results on the logarithmically enhanced QED 
corrections to the forward backward asymmetry and to the decay width in the high-g^ region 
and also discuss some nonperturbative subtleties. Furthermore, we analyse the relation of 
the collinear logarithms and the experimental angular cuts. In Section H] we present the SM 
predictions for the branching ratio and forward-backward asymmetry. In Section [5] we draw 
our conclusions. 



2 Master formulae 

In this section, we present master formulae for the various observables. We closely follow here 
the notations of Ref. [30] in order to keep the paper compact. The new results are explicitly 
derived in the next section. 



2.1 Dilepton invariant mass spectrum 

Following Ref. [30] , the master formula for the B Xsi~^i~ branching ratio reads 



dB{B ^ Xs 



ds 



8{B X,eu] 



exp 



cb 



4 ^aCs) 



C 



where s = q'^/mb^pde- and are defined by 



dT{B Xs 



ds 

T{B X^eu) 



(3) 



(4) 
(5) 



and we have = 1 + (9(q;s, Oem, A^/m^). The normahzation with the measured semileptonic 
decay rate minimizes the uncertainty due to the fifth power of m^; the factor C, 



C 



uh 



ch 



(6) 



is used to avoid spurious uncertainties due to the h X^ev phase-space factor. This factor 
can be determined from a global analysis of the semileptonic data [31]. The quantity 
can be expressed in terms of the low-scale Wilson coefficients and various functions of s that 
arise from the matrix elements. The main formula reads 



(7) 



where C^^{fib) 7^ Ci{fib) only for i = 7,8. The functions s) are given in Eq. (116) of 

Ref. [30] and can be expressed in terms of the coefficients listed in Table 6 of Ref. [30], as 
well as the building blocks Sgg, S'77, 5*79, and Sioio given in Eqs. (112-115) of Ref. [30]. 

In the low-g^ region, there are only a few modifications compared to the analysis in Ref. [30]: 
we update the parametric inputs, the influence of the cc resonances are taken into account and 
the two-loop matrix element of the operator Pg is given in a more precise way (see Eqs. ( HHi) 
and f H9l) of the present paper). The combined effect of these improvements does not change 
appreciably the central value for the branching ratio obtained in Ref. [30]. 

In the high-g^ region, there are several additional changes necessary besides the ones men- 
tioned above. The electromagnetic contributions to the building blocks Sqq, Sjt, Sjq, and 5'ioio 
are different and are presented here for the first time in Eqs. pUj) - fH7j) . 0(l/m^) corrections 
are numerically relevant and explicitly given in Eqs. (6) of Ref. [24]. Finally, the functions F^, 
F2, Fj, F2, and F^ that appear in the coefficients (listed in Table 6 of Ref. [30]) are 
different in the high-g^ region and are known only numerically. They have been calculated in 
Ref. [11] Q 



2.1.1 New normalization 

The authors of Ref. [24] have shown that it is possible to drastically reduce the size of 1/ml and 
l/ml power corrections to the integrated decay width, by normalizing it to the semileptonic 
b —>■ uiu rate integrated over the same g^-interval. This procedure will help reducing the 
uncertainties induced by the large power corrections to the decay width integrated over the 
high-s region. The new observable is defined as follows: 



7^(so) 



1 ^^dT(B ^ XJ+e-) 



so ds 



1 ,.dr(50 ^ 



so as 



XJu) 



VtsVtb 



ub 



2 fl 



4 ^ii{s) 



(8) 



^Numerical expressions of these functions can be obtained from the authors upon request. 
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where ^u{s) was defined in Eqs. and (jl]) and the differential is given by 



dr, Gj,\Vub\m: 



6, pole 



dS 192.3 I!*) 

Exphcit expressions for the 0(1, a^, contributions to can be found in 

Eqs. (5,11) of Ref. [24]. The O(a^) correction to ^u{s) is given by q;^X2(s) with X2 defined in 
Eq. (60) of Ref. [32]. We have not included electromagnetic corrections of order 0{k,) in the 
normalization because they are presently not known. However, the impact of those effects in 
the conventional approach where the fully integrated semileptonic decay rate is used as nor- 
mahzation is less than 2.5%. We would also like to note that in Eq. ([H]), contrary to Eq. (jS]), 
corrections of 0{asK,) to $u(s) can contain residual logarithmically enhanced terms once the 
integration over s is restricted to the high-g^-region. We leave the investigation of terms of 
0{k,) and 0{asn) to a future analysis. 

Note that in the definition of TZ we used only neutral semileptonic decays in order to reduce 
the uncertainties coming from weak annihilation contributions (see further remarks in Section 
4.1). 

Obviously, for the comparison of the theoretical ratio TZ with experiment, a separate mea- 
surement of — * Xuii^ and X^lv with a high-g^ cut is necessary. These are quantities 
with larger experimental uncertainties. However, they are not expected to negate the drastic 
reduction of the theoretical uncertainty in the ratio TZ (see section W?I\ for numerical results). 



2.2 Forward-backward asymmetry 



In complete analogy to the formula for the branching ratio, one can derive also a formula that 
expresses the FBA (see [34]) in terms of the low-scale Wilson coefficients and various building 
blocks. Normalizing to the semileptonic B X^ev decay width, we get 



AAfb{B -> X, 



where $^^"^(5) is defined by {z = cosOi 



B{B X,eu] 



cxp 



V*Vt, ' 4 <|.f/^( 



c 



dz sgn{z] 



d'^T{B X, 
dsd^; 



-1 



5 



(10) 



^FBA 



Re[cf (^,)Cf*(/x,)ff,,(/i5,s) 

i<3 



12) 



Again, the functions Hij{fii,, s) depend on the coefficients Mf^ listed in Table 6 of Ref. [30] and 
on the building blocks 5*710 and S'910: 

E Re(M^M/o*) 5aio + A/fii , when z = j 

A=7,9 

(13) 



E (M^Mf* + Ml^Mf*) Saio + ^Hij , when i j 
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where S'710 and Sgio are given by 



'710 



jgio 



-6(l-s)2{l + 8 5s [f7io{s)+u 



,(1) 



(em) 
710 ' 



s) +U 



(em) 



/710(s' 



-16 5.2 



n(2) + 4M«/no(S)]} 



A, 



8 Ai Ao , 

^ s- 6 ^ 1 - 14s + 9s^) 

-3s(l-s)2{l + 8«s [hw{s)+u 



(14) 



,(1) 



+ KU 



(em) 



+ 8 ttsK 



+ lQal [m(2) + 4m«/9io(s)]} 



4Ai 



mt 



12 Ao o / N 
— ^ 4 - 3s . 
mi 



(15) 



The functions /jio(s) (^ = 7,9) can be found in Eqs. (15-17) of Ref. [9]. The functions cufj™'' 
represent the electromagnetic corrections to the matrix elements which are calculated for the 
first time in the present paper (see Eqs. ([H2H3H1))- The Sab also include non-perturbative 
0(1/7715) corrections from Refs [19,22]. Contrary to the expression for the branching ratio the 
quantity Ai, which is related to the kinetic energy of the 6-quark, does not drop out here. The 
quantities 



AK. 



(16) 

have the same meaning as in Eq. (117) of Ref. [30]. They need to be included only for i = 1,2. 
The additional In (777^ /?7i|) -enhanced electromagnetic corrections for the FBA read 



6210 = ~2AHl--sray4\^\s) 
_ 4 

Clio — 06210- 



(17) 



The function ci;2i™'*('5) is again new (see Eq. (IMl) ). while the 0{AQQY)/ml) non-perturbative 
contributions were calculated in Ref. [21] 



C210 

Clio 



A2 



3ml 



:i-s)2(l + 3s)F(r) 



~fi C210 , 



1^ 



where r = I/t/c = ■s/(47ri^) and the function F{r) is listed for example in appendix A of 
Ref. [30]. The finite bremsstrahlung contributions bij were calculated in Ref. [13]. We do not 
present these corrections here but do include them in the numerical analysis. 



2.3 Normalization of the FBA 

As discussed in section 2, the normalization of the FBA by the semileptonic or by the B —>■ 
Xsi~^i~ decay rate is important in order to cancel the overall factor ml^^^^. Moreover, the 
elimination of renormalon ambiguities requires the analytical conversion to a short- distance 
mass and a complete expansion in ag [35,36]. 
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If we extract the zero of the FBA without using any normahzation, the conversion and the 
expansion of the overall ml ^^^^ factor induces a large dependence on the scheme we choose for 
the short distance b mass. In fact, we find a — ll%(+5%) shift in switching from the IS scheme 
to the MS (pole) one. This large scheme rrib dependence was overlooked in previous analyses of 
the zero. It can be eliminated by normalizing the FBA by a rate proportional to ml and fully 
expanding in a^. 

The most natural choice is to normalize the FBA directly by the differential B X^i^i^ 
differential decay rate because this normalisation is also used in the experiment. Unfortunately, 
this option is not viable because of the poor convergence of the perturbative series for the latter 
(the most serious problem is that the lowest term in the expansion is sub-dominant). It is 
therefore not surprising that also this procedure leads to a large mj, scheme dependence of the 
zero: we find a shift of —12% when one compares the IS to the MS scheme. 

Therefore, we propose here a definition of the FBA that, at the same time, removes the 
ml dependence, does not suffer of renormalon ambiguities, has a nicely behaved perturbative 
expansion, and can directly compared with the experimental quantity. The following double 
ratio satisfies all the requirements: 



where the quantities are defined in Eqs. ([3]) and (ITOl) . The critical point is that each of the 
ratios in brackets gets fully expanded in as, but no overall expansion is done. This is allowed for 
two reasons: on the one side, the fully expanded numerator and denominator are renormalon 
free, on the other one, they separately correspond to two measurable quantities. As we show 
in Sec. 14.51 also the scheme dependence of the zero turns out to be minimal (the difference 
between the IS, MS and pole schemes is about 1%). 

This discussion shows that the estimation of the perturbative error of the zero via the com- 
parison of different normalisations as done in [12] is questionable because some normalizations 
are affected by renormalon ambiguities and the problem of the large scheme dependence of the 
b quark mass should be addressed. As we will discuss in section 14.51 the standard analysis of 
the scale dependence offers a better guidance to the understanding of the perturbative errors 
in our predictions. 

2.4 Perturbative expansion 

Some remarks about the perturbative expansion are in order: Large logarithms of the form 
ttg log(mfe/Mvi/) have to be resummed at a given order in perturbation theory using renor- 
malization group techniques. In the case of 6 — si^i~ , the first large logarithm of the form 
\og{mb/Mw) arises already without gluons. Moreover, the amplitude of 6 — *• si^i~ is propor- 
tional to aem- So naturally we have an expansion in as and in k = aem/c^s while resumming all 
powers of aslog{mb/Mw)- The Leading Order (LO) are of order k,, NLO are proportional K,as, 
NNLO are proportional na'^. 

It is well-known that this naive expansion is problematic, since the formally-leading 
0(1/ as) term in Cg is accidentally small and much closer in size to an 0(1) term. In addition 



ds ds 




(s) <I>a(s) 



(19) 
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the NLO terms are enhanced by a factor of mf / {M^^r sin^ 9w) ■ Therefore, also specific higher 
order terms in the general expansion are numerically important. 

The b si~^i~ decay amplitude has the following structure (up to an overall factor of Gp)'- 



A 



A 



LO 



as Anlo + c?s Annlo + 0{p?s, 

A em I ^ /I em i ^ 2 /i em , /n ^ 3 ^ 



(20) 



with Alo 



a. 



A 



NLO 



and A 



LO 



a. 



N LO- 



AM these terms are included in the numerical 



analysis in a complete manner, together with the appropriate bremsstrahlung corrections, while 
also the term Annlo is practically complete due to the calculations in Refs. [5-14]. The 
only missing parts originate from the unknown two-loop matrix elements of the QCD-penguin 
operators whose Wilson coefficients are very small. 

Among the contributions to A'^^nlo^ include only the terms which are either enhanced 
by an additional factor of mf / {M^^ siia^ 9w) (with respect to AJIJI^q) [12] or contribute to the 
ln(m^/m£)-enhanced terms at the decay width level. As discussed above, the latter terms 
were calculated in Ref. [30] for the dilepton mass spectrum in the low-g^ region, while the 
corresponding terms for the high-g^ region and the FBA are calculated for the first time in the 
present paper. 

The perturbative expansion of the ratio $£^(s)/$u has a similar structure to that of the 
squared amplitude (up to bremsstrahlung corrections and nonperturbative corrections) 



A' 



A'lo + '^s '^AloAnlo + Oil {A%i^o + '^AloAnnlo) 
+al 2{AnloAnnlo + •.•) + 



em 

NLO-^LO 



2ALoAro + as2{A^r.nA 

2{AnloA'}^lo + A 
+al2{ANLoA7rNLn + A 



^LO-^NLO) 



A em 

NNLO-^LO 



+ A,^oA\ 

+ Ai^qA^^nlo) 

AT,.n + ...)+ (^(«') 



em 

NNLO^NLO 



(21) 



We assume that all products in Eq. ([7j) are expanded in and k. Regarding QCD, a strict 
NNLO calculation of $«(s)/$„ should only include terms up to order K^a^. 

In the numerical calculation of $££(s)/$u, however, we include all the terms that are written 
explicitly in Eq. (1211) . The term Anlo Annlo of order K^a^ is formally a NNNLO term, 
but numerically important. Within the electromagnetic corrections the same is true for the 
term Af^p^oANNLO- We emphasize that those terms beyond the formal NNLO level which are 
proportional to jCyp and jCgp are scheme independent. The dots in Eq. fl?Il) stand for the 
unknown terms AloAnnnlo 



and AJ^A 



mNNLO 



and, consequently, can safely be neglected due 



to A 



LO 



asA_ 



NLO: 



A 



LO 



and OsAnnnlo 



< A 



NNLO- 



Thus, one can argue that our 
calculations pick up the dominant NNNLO QCD corrections. In the following we will therefore 
call the accuracy of our calculations improved NNLO. § 



^In fact, this scheme is very similar to the one proposed in Ref. [6] and used in many previous analyses of 
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Finally let us comment on the mass schemes that we use: The pole mass of the b quark 



hindered by renormalon ambiguities that appear in the relation between any short distance 
mass definition (e.g. MS, IS, kinetic schemes) and the pole one. We eliminate these renormalon 
uncertainties utilizing the Upsilon expansion described in Refs. [35,36]. Every occurrence of 
the pole mass is converted analytically to the IS*- mass before any numerical evaluation of the 
branching ratio is performed. In our analysis we use the conversion formula up to order [36]. 
We follow a similar approach for the treatment of the charm pole mass that appears in the 
calculation of some matrix elements. We adopt the MS charm mass as input and expand the 
pole mass at order using the formulae presented in Ref . [36] . For what concerns the top mass, 
we take the pole mass as input and convert it to the MS scheme at order a^. The perturbative 
expansion of the FBA proceeds along the same lines. 



3 Analytical results 

In this section, we calculate the electromagnetic corrections to the matrix elements for the FBA 
in the low- and high-g^ region and for the dilepton mass spectrum in the high-g^ region and 
discuss the phenomenological implications of those coUinear logarithms in the ee and fifi final 
states. 

Moreover, we present a more precise fixing of the two-loop matrix elements of the Pg operator. 
Finally, we make some remarks on some non-perturbative subtleties. 

3.1 Log-enhanced corrections to the FBA, low- and high-g^ region 

First we derive the basic formulae to obtain the expressions for the ln(m^/m|)-enhanced correc- 
tions to the FBA, where we shall focus here on the un-normalized asymmetrj{§ corresponding 
to the LHS of Eq. (ITTI) . We follow closely the notation of Section 5 of Ref. [30]. 

It was shown in the Appendix of Ref. [19] that the angular FBA with respect to O^, the 
angle in the dilepton c.m.s. between the directions of the momenta of the decaying B and the 
positively charged lepton, is equivalent to the energy asymmetry between the two leptons in the 
rest-frame of the decaying B. Events in which cos 6'^ > in the dilepton c.m.s. correspond to 
events in which E_ > measured in the 5-meson restframe, where E± denotes the £^-energy. 



the b s£^£~ observables. In these works the formaUy-leading, but accidentally small 0{l/as) term in Cg is 
treated as 0(1) and absorbed into the NLO coefficient: as consequence the two-loop matrix element of Pg and 
the three-loop mixing of the four-quark operators into Pg are left out even though they are formally NNLO 
terms. This is the main difference between the NNLO scheme proposed in Ref. [6] and the improved NNLO 
used in the present manuscript. 

§In this subsection we suppress the normalization factor, thus, the symbol Afb denotes the un-normalized 
FBA, see Eqs. ([23l) . ([27| . and ([3T|) . 
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Sticking to the latter frame and defining the scaled energies 

y± = ^. 22 

we can write the fully differential FBA as 

dApB = PF ds dy+ dy^ 5{l + s - y+- y_) \Af sgn(y_ - ?/+), (23) 
with the pre-factor 

327r3 • 

We are mainly interested in ln(mf/m|)-enhanced electromagnetic corrections. They are derived 
by means of the splitting function and we shall adopt the kinematics from Figure 2 of Ref. [30]. 
In the collinear limit the fully differential FBA reads 

dAplcon = PF dx ds dy+ dy^ 6{1 + s - y+ - y^) f^;^\x) \Af sgn(|/_ - y+) . (25) 

We shall only retain the In (m^/m|) -enhanced part of f^"^\ which then becomes independent of 
E, 

/WW=4a.Il±(l-^£lluf!^^) . (26) 

The approximation of dropping the non-log-enhanced terms in the splitting function is justified 
by the fact that even for high the logarithmically enhanced part dominates numerically. This 
can be seen, for instance, by the fact that the non-log-enhanced term of (1 — s)^ (1 + 2s) u!^(^\s) 
(see Eq. (94) of Ref. [30]) vanishes at s = 1. 

In the squared amplitude in Eq. ( l25l) we keep only those terms that are relevant for the FBA, 
i. e. all terms that do not drop out upon integration over the sign-function, 

1^1' = 2Re[C^oC; + a,l^C,o{C2 + CFC,)f;{s)]{P^o)tree{P9)Le 

+2Re[C,oC;] (Pio)trcc(P7);ee • (27) 

As in the case of the decay width we must consider the difference 



"■^Fg,coll,2 _ "-^FB,coll,3 

ds ds 



{2i 



where we stay differential in the double and in the triple invariant, respectively. At this point 
more care is required compared to the calculation of the decay width. Due to the emergence 
of the sign-function in Eq. fl23|) we must distinguish between photon-emission from the ^+ and 
from the l~ in case of the double invariant. In the former case we have to stay differential in 
s = {xpi + p2y/ml, in the latter in s = {pi + xp2)'^ /ml. Also the y± change accordingly. We 
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therefore have 



dA 



(m) 

FB,coll.2/± 



ds 



±PF 



— J^T J 

1— VS x(l + x) 

I 

S/x 



1-1 1 



1 



J dx J dy± + J dx jdy± — j dx j dy± 

s/x 1-v/I s/x 



X 



x] 



X 



*s/x ; yif:^l-y±+s/x 



(29) 



The two expressions corresponding to upper and lower sign should be equal due to the antisym- 



metry of \Af in ^ y_. The case of the triple invariant is simpler since we stay differential 
in s = (pi +P2)^/^l, namely 



dA^""^ 
ds 



±PF dx 







1+J 

2 



Jdy±- J dy± 



1+1 

2 



y^^l-y±+s 



(30) 



where the expression for the upper and lower sign should again be equal due to the antisymmetry 
of \Af. We finally have to combine the expressions according to Eq. 
The corrections to the unnormalized FBA read 

dAApB 



ds 



^M\v;xf A$f/^(s) 



with 



(em) , 
710 y 



In 



487r3 
-24 s 

m3 



VttV;f{l-sfaM -48 



(cm) , 



Re [^9^*0] 4i?(5) + «s« Re \{C2 + CpC^)Cl, 4\'^\s] 
+ ln(l- 



(31) 



7 - 16 Vs + 9 s 



LU, 



(em) I 
910 I 



4 (1 - s) 

5 -16 73 + lis 
4 fl-s) 



s Ins 



l-s) 



(32) 



l-5s ^^/ 1 + Vs 



1 - s 



+ ln(l - VI) 

(1- 3s) Ins 
(l-s) 



^210 \S} 



In 



S7(s) +2S((s) 
24 5(1-5)2 



(33) 



(34) 



The functions uj^^\s) and C(j9'[™''(s) are known in the entire g^-region, whereas the function 
^2\o\^) 'was obtained by a least squares fit. The function S7, valid in the low-s-region, reads 

S7(s) = -0.259023 - 28.424 s + 205.533 5^- 603.219 5^ + 722.031 sS (35) 
T.\{s) = [-12.20658 - 215.8208(5 -a) +412.1207 (s-a)^] (s-a)2^(s -a) , (36) 
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with a = {Aml/miy . In the high-s-region the function Sy reads {6 = 1 — s) 

£7(5) = 77.0256(52- 264.705 5^ + 595.814 5^ - 610.16375^ 
= 135.858 (5^- 618.990 5^ + 1325.040 (5^- 1277.170 (5^ 

The polynomials in the high-s-region were obtained such as to have a double zero at s = 1 



(37) 
(38) 



3.2 Log-enhanced corrections to the BR, high-g^ region 

The log-enhanced corrections to the differential decay width given in Eq. (jlj) read 



^ ''■\VtbV;fil-s)^aM 8(1 + 2:3) 



\^io\ ^1010 1-5 



ds 487r3 

+a,K Re [(^2 + CfCi)C; 47^^)1 + 5,^^ (C2 + CpCiY 47^(5) 



+96 



a,K Re [CjC;] ujif'\s) + a'tn' Re (C2 + CfCi)C; cj^r^(s 



y* (cm) , 



+8(4 + pc5,V|C7p4r^(S) 



(39) 



Exact analytical expressions are known for the functions ujg'(^\s), co'ioio(s), u!^^-^\s), and 
4r^(s) (see Eqs. (94) and (100) - (102) of Ref. [30]). Therefore, they hold in the entire 
g^-region, while the other w-functions are obtained by a least-squares fit in the high-g^-region 
(for fixed values of rUb and rric): 



(em) , 



29 



(em), 
27 I 



= In 


/m| 




\m\ 






= In 


/m| 




Vm| 




64 


+ 




81 






= In 


/m| 







YjAs) + i Si(s 



16 



+ ^^S(^) In 



5GeVy 



8(l-s)2(l + 2s) 

8(l-s)2(l + 2s) ' 9(1 - s)2(l + 2s) ""V5GeVy 



(40) 



In 



(em) , 
1010 I 



V5GeV 



9 



96(1 - s)2 

with the functions Sj, valid in the high-s-region ((5 = 1 — s): 

-148.061 (5^ + 492.539 - 1163.847 5^ + 1189.528 , 
-261.287 52 + 1170.856 5^ - 2546.948 (5^ + 2540.023 b^ 



S4(s) 
Eifs) 



S5(5) 

S6(s) 



-221.904 (5^ + 900.822 b^ - 2031.620 5^ + 1984.303 b^ , 
-298.730 (5^ + 828.0675 b^ - 2217.6355 b'^ + 2241.792 b^ 
-528.759 (5^ + 2095.723 5^ - 4681.843 b'^ + 5036.677 (5^ . 



(41) 
(42) 



(43) 
(44) 
(45) 
(46) 
(47) 
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The fits are excellent for s > 0.65. Again, the polynomials were obtained such as to have a 
double zero at s = 1. 

3.3 Phenomenological implications of the collinear logarithms 

After inclusion of the NLO QED matrix elements calculated in the previous subsections the 
electron and muon channels receive different contributions due to terms involving ln(m^/m|). 
We emphasize that this is the only source of the difference between the electron and muon 
channel. 

However, as already pointed out in Ref. [30], the presence of this logarithm is strictly related 
to the definition of the dilepton invariant mass. The collinear logarithm ln(mf/mf) disappears 
if all photons emitted by the final state on-shell leptons are included in the definition of g^: 
{Pii + Pi2Y ~^ ipii + Ph +^7)^- If none of these photons was included in the definition of - 
i. e. if a perfect separation of leptons and collinear photons was possible experimentally - then 
our expressions containing \n{ml/mj) were directly applicable. This is the case for muons since 
the separation of muons and collinear photons is practically perfect [46] . For electrons with the 
current Babar and Belle setups the lepton mass gets replaced by an effective mass parameter 
A which is found to be of the same order as m^. Hence our numerical results in section H] for 
muonic final states should be applied also to the case of electronic final states in the present 
Babar and Belle setups. All results in section H] for the electron channel are given under the 
assumptions of perfect separation of electrons and collinear photons. We refer the reader to 
section 6 of Ref. [30] for more details. 

3.4 Two-loop matrixelement of Pg 

In Ref. [12], an estimate of the two-loop (NNLO) matrix element of the operator Pg - denoted 
ujgg {s) in Eq. (74) of Ref. [30] - was presented which was neglected in all previous analyses. The 
estimate is based on the fact that the QCD corrections to 6 — si'^i' are identical to those of 
b uiu {t — i> Mu) , in the limit of vanishing strange (bottom) quark mass. However, the specific 
two-loop calculation of the b uiu decay was available as an expansion in (1 — s) [37], while 
the two-loop contribution of the top quark decay was only known as an expansion of M^/m^ 
(which translates into an expansion in q^/ml for b —>■ sit) up to the second order based on 
Fade approximation methods [38]. In the meanwhile, the QCD corrections to 6 ^ ulv are also 
known as an expansion in s [32,39]. 

(2) 

Thus, in the low-s region we approximate the function cugg by the m s = u expanded results 

= X2{uj = s)/X^{u = s) (48) 

where X2{uj) is given in Eq. (60) of Ref. [32] and Xq{uj) in Eq. (2) of Ref [39]. 

Analogously, in the high-s region we approximate the function cjgg by the in (1 — s) = 5 
expanded results of ^2(5 = 1 — s) and Xq{5 = 1 — s) as given in Eq. (2) and (3) of Ref. [37]. 
Note that the normalization of the Xj differs by a factor of 2 in Ref. [37] and Refs. [32,39]0 

^There is a typo in Eq. (4) of Ref. [37]: the powers of the last four (5's have to be increased by 1 each. 
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3.5 Some nonperturbative sublet ies 

For a detailed discussion of nonperturbative corrections we refer the reader to Section 5 of 
Ref. [11]. Here we add few remarks on some specific issues. 

The long-distance corrections related to the cc intermediate states originate from the non- 
perturbative interactions of the cc pair in the process B — > XgCC — > Xsi~^i~ ■ If the dilepton 
invariant mass is near the first two J^'-'' = 1 cc resonances and ^'), this effect is very 
large and shows up as a peak in the dilepton mass spectrum which can easily be eliminated by 
suitable kinematical cuts. More delicate is the estimate of such long-distance effects away from 
the resonance peaks in the perturbative windows. 

Within the KS-approach [25], one reabsorbs charm rescattering effects into the matrix 
element of Pg; the effects of b ^ ccs operators is estimated by means of experimen- 
tal data on (T(e+e^ — > cc hadrons) using a dispersion relation. To be more specific, the 
function h{z,s) appearing in Eq. (A. 5) of Ref. [11] (which corresponds to the expression 
diVc) + 8/91og(mb/mc) - 4/9 in Eq. (72) of Ref. [30]) is replaced by {z = ml/ ml) 

h{z, s) hiz, 0) + -P / ^ ds'-f^ + z-K^{s) , (49) 

where R^^^{s) = (j{e^e~ cc) / a{e^e~ — > ii^ix~). This method is exact only in the limit in 
which the B X^cc transition can be factorized into the product of sh and cc colour-singlet 
currents. It is possible to take into account non-factorizablc effects and reproduce the correct 
hadronic branching ratios by multiplying P^ad by a purely phcnomcnological factor n = 2.6. 
However, a model-independent way to estimate the non-factorizable cc long-distance effects far 
from the resonance region exists by means of an expansion in inverse powers of the charm- 
quark mass [20,21]. Having included those l/rric corrections, an 'inflation' of the factorizable 
KS-corrections using the phenomenological factor k would lead to a double- counting. The 
numerical impact of the KS contribution to the integrated rate in the low- and high-s regions 
is about -1-2% and —10%, respectively. 

In Section 5 of Ref. [11], another large nonperturbative error for the high-g^ region was 
located in the linear l/rniy corrections. The physical observables defined in terms of a g^-cut 
are sensitive to a l/m^ term via the the uncertainty on the value of m(,^poie or equivalently via 
the relation of the pole mass and the hadron mass 

MB^m^{l+K/mj, + 0{l/ml)). (50) 

Choosing as reference cut the value Smin = 0.6, the integrated normalized dilepton mass spec- 
trum defined in terms of g^j^ can be written as 

X f~^f~) 

-n-cutWminJ 




= n -6.2(^-0.6^ +0 f^-0.6^ \xf\sR{s), (51) 



0.6 
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as{M,) = 0.1189 ± 0.0010 [40] 


me = 0.51099892 MeV 


aeiM,) = 1/127.918 


= 105.658369 MeV 


= sinHw = 0.2312 


= 1.77699 GeV 


\VtsVttW = 0.962 ±0.002 [41] 


me(mc) = (1.224 ±0.017 ±0.054) GeV [42] 


\VtsVtb/Vub\^ = (1.28 ±0.12) X 10^ [41] 


ml^ = (4.68 ± 0.03) GeV [31] 


BR{B Xcez/)exp = 0.1061 ± 0.0017 [43] 


mi,poie = (170.9 ±1.8) GeV [44] 


Mz = 91.1876 GeV 


TUB = 5.2794 GeV 


Mw = 80.426 GeV 


C = 0.58 ±0.01 [31] 


Af = (0.12 ±0.02) GeV^ 


Pi = (0.06 ± 0.06) GeV^ [31] 


Af = (-0.243 ± 0.055) GeV^ [42] 


/o±/, = (0±0.2) GeV^ [24] 


/°-/, = (0±0.04) GeV^ [24] 


ft = (0 ± 0.4) GeV^ [24] 



Table 1: Numerical inputs that we use in the phenomenological analysis. Unless explicitly 
specified, they are taken from PDG 2004 [45]. 



which implies 



6R 



cut 



7.4 



Snib 



(52) 



-Rcut mb 

Using the pole mass scheme with Srub = 0.1 GeV, this leads to a ~ 15% error on -Rcut- However, 
this error gets now significantly reduced in our updated analysis using the kinematical IS scheme 
for the nib mass. 



4 Numerical results 

The numerical inputs that we adopt are summarized in Table [TJ 



4.1 Branching ratio in the high-g^ region 

For the branching ratio integrated over the region > 14.4 GeV^ we find: 



Kjhigh 

MM 



2.40 X 10"' 1 ± 



+0.01 
-0.02 



± 



MO 



+0.14 
-0.06 



± 0.02^, ± 



Mi) 



+0.006 
-0.003 



C.nic 



± 0.05,„, ± 



+0.0002 
-0.001 



±0.002cKM ± 0.02br,, ± 0.05a2 ± 0.19^, ± 0.14/, ± 0.02/„ 
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2.40 X 10-^ > (53) 

±0.05. 



B^'J^'- = 2.09 X 10-M 1 + ^^4 + ^0 08 ±0.02^,+ 

J 111, 



+0.005 
-0.0009 



+0.0003 
-0.002 



Mo 

±0.002cKM ± 0.02br,, ± 0.05a, ± 0.22,, ± 0.16^^ ± 0.02^,^^ 
= 2.09x10-^(11^2). (54) 

The scale uncertainty has been estimated by varying the matching scale /xq and the low-energy 
scale fib by factors of 2 around their central values (120 and 5 GeV, respectively). In our 
approach we normalized the b sii decay width to the semileptonic rate averaged over neutral 
and charged B mesons; hence, we average and f^. From the values quoted in Tabled] we 
obtain: = ± 0.2 and = ± 0.1 1 

The other parametric uncertainties are obtained by varying the inputs within the errors 
given in Table [H We assume the errors on C and rric to be fully correlated. The total error 
is obtained by adding the individual uncertainties in quadrature. We note that here and in 
the following all errors are parametric or perturbative uncertainties only and that subleading 
nonperturbative corrections of order 0{asA/mb) will give an additional uncertainty. 

Log-enhanced QED bremsstrahlung corrections shift the central value by about —8% and 
—20% for the muonic and the electronic final state, respectively. These relative shifts are much 
larger than the corresponding ones in the low-g^ region. In fact, these corrections vanish when 
integrated over the whole dilepton invariant mass spectrum, but are relevant if one is restricted 
to certain regions in phase space. One observes that the absolute shifts in the central values 
of the branching ratio are of similar size but of opposite sign in the low-g^ compared to the 
high-g2 region, namely for the muonic final state +3 x 10"^ and —2 x 10"^ respectively. But 
since the differential decay width decreases steeply at large g^, the relative effect is much more 
pronounced in the high-g^ compared to the low-g^ region. See also Fig. [T] for illustration. 

4.2 Ratio 7^(so) 

The numerical results we obtain for the ratio TZ{sq) discussed in Section [2.1.11 are for sq = 14.4 
GeV^: 

^(so)Jf = 2.29 X 10-3 (l ± 0.04scaic ± 0.02™, ± 0.01c,™. ± 0.006^, ± 0.005,, ± 0.09ckm 
±0.003a, ± 0.05p, ± 0.03^0+^, ± 0.05^o_;, 



"Regarding the weak annihilation contributions (/°, /+, , /+) we follow the analysis in Ref. [24]. Here are 
some additional remarks in order: Using flavor SU{3) one concludes that /° ~ However, we stress that the 
additional assumption /" ^ /+ comes from neglecting iso-singlet effects in the SU{3) decomposition of bs sb. 
Nevertheless, due to a mismatch in the flavor of the spectator quark in the B meson with the quark appearing 
in the operators ( bu ub and 6s sb), it is not unreasonable to assume f^,fu<fu- This leads to the relation 
/° ~ /° ~ /+ < /+ [24] that we use in the numerics. We note again that in the definition of TZ we use only 
neutral semileptonic decays, thus, the dependence on disappears in the ratio. The actual numerical inputs 
we adopt are taken from Ref. [24,33]. Because all our expressions for the branching ratio and the ratio TZa are 
linear in these non-perturbative parameters, hence it is trivial to adjust our results to accommodate for different 
sets of non-perturbative inputs. 
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Figure 1: The log-enhanced part of (1 — s)^(l + 2s) ujg(^\s), see Eq. (94) of Ref. [30]. This 
function integrates to zero, but gives contributions of opposite sign to the low-s and high-s 
branching ratio, respectively (shaded areas). The moduli of the absolute shifts are of comparable 
size. However, the relative effect is much larger in the high-s compared to the low-s region due 
to the steep decrease of the differential decay width at large s. 

= 2.29 X 10-^(1 ±0.13) , (55) 

^(So)e?'' = 1-94 X 10-^^(1 ± O.OGscale ± 0.02^, ± 0.02c,rr.c ± 0-004^, ± 0.006a, ± 0.09cKM 

±0.01,, ± 0.09,, ± 0.05;o+;, ± 0.05;o_^,) 
= 1.94 X 10^^(1 ±0.16) . (56) 

Note that uncertainties from poorly known 0(l/m^) power corrections are now under control; 
the largest source of error is Vub- 

4.3 Branching ratio in the low-g^ region 

For the branching ratio integrated in the range 1 GeV^ < m^^ < 6 GeV^ it was found in Ref. [30] : 



1.59 ± O.OSscaie ± 0.06„, ± 0.024c,^^ ± 0.015^, 

±0.02«,(M^) ± 0.015cKM ± 0.026br,,] x 10"^ = (1.59 ± 0.11) x 10"^ , (57) 
1.64 ± O.OSscaie ± 0.06^, ± 0.025c,m. ± 0.015^, 
±0.02, 



.^^^M,) ± O.OlScKM ± 0.026brJ X lO-'' = (1.64 ± 0.11) x 10"^ . (58) 

We note that we find for B^^ a +1.8% shift of the central value due to the KS-corrections which 
were discussed in Section 3.4 and which were not included in the previous analysis. Moreover, 
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the update of the input parameters (CKM, mt, as) leads to a —3.1% shift; thus, we end up to an 
overall change of —1.3% within our new phenomenolgical analysis compared with the previous 
one in Ref. [30]. The total error, being from parametric and perturbative uncertainties only, 
remains unchanged. 



4.4 Integrated FBA in low-s region 

By the end of the current B factories the fully differential FBA will not be accessible experi- 
mentally, contrary to the integral over one or more bins in the low-g^ region. However, these 
integrals can already serve of gain important information on the shape of the FBA and to 
constrain the parameter space of new physics models (see for instance [47]). We now give the 
results for the integrated FBA based on Eq. ffT^ . We subdivide the low-s region into the two 
bins s e [1,3.5] GeV^ and s E [3.5,6] GeV^, which we will call bin 1 and bin 2, respectivel\l**1. 
We then integrate the numerator of the double ratio, Eq. ( fTOl) . over the respective bin, and 
divide the result by the denominator of the double ratio, Eq. ([3]) integrated over the same bin, 
i. e. we first integrate FBA and branching ratio separately and subsequently divide the two 
numbers. Our results read: 



(•^^?)bi 



binl 



binl 



bin2 



/bin2 



-9 


09 ± 0.83scalc ± 


03m,i 


± 0.24„,,c 


± 0.20™, 


± O.lSa^Mz) 


± 0.02a2 


-9 


09 ±0.91 


%, 












-8 


14 ± 0.80scale ± 


03mj 


± 0.23^^,0 


±0.19™, 


± 0.18a,{Mz) 


± 0.02a, 


-8 


14 ±0.87 


%, 












+ 7 


80 ± 0.55scale ± 




±0.31^,,c 


± 0.34™, 


± 0.16a4A/z) 


± 0.20a/ 


+ 7 


80 ±0.76' 


%, 












+ 8 


27 ± 0.47scaic ± 


02 


± 0.30^,,c 


± 0.33™, 


± 0.15a,{Mz) 


±0.19a2 


+ 8 


27 ±0.69' 


% . 













(59) 

% 

(60) 

% 

(61) 

% 

(62) 



The total error is obtained by adding the individual ones in quadrature. For the entire low-s 
region, we get 



/ low 



% 



low 



- 1.50 ± 0.783caic ± 0.02™, ± 0.29™,,c ± 0.27™, ± O.IS^jm,) ± 0.10a, ^ 

- 1.50±0.90] % , (63) 

- 0.86 ± 0.73scaic ± 0.01™, ± 0.28™,,c ± 0.26„, ± O.IS^^^Mz) ± 0.10a, ] % 

- 0.86 ± 0.85] % . (64) 



The relative errors in the respective bins are considerably smaller than for the entire low-s 
region since the respective values in each bin are similar in size and of opposite sign. 



"Predictions for different bins can be produced upon request. 
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IS 


MS 


pole 


/i 


3.50 


3.47 


3.52 


e 


3.38 


3.34 


3.41 



Table 2: Dependence of the zero of the FBA on the b quark mass scheme. The input values 
are mb{mb) = 4.205 GeV [48] and m6,poie = 4.8 GeV. 



4.5 Analysis of the zero of the FBA 

The basic formula for the extraction of the zero is Eq. fll9p . We expand the two ratios separately 
in as, n, Ai and A2, and keep all the terms as specified in section [231 It is understood that also 
the conversion of the mass scheme for the botton, the charm and the top quark is performed 
in the way described there. The results for q^, the zero of the FBA in the low-s region, are 



3.50 ± O.lOscaic ± 0.002„, ± 0.04^^,c 

±0.05„, ± 0.03«4A/z) ± O.OOIai ± O.OIa^ 1 GeV^ = (3.50 ± 0.12) GeV^ 



(65) 



3.38 ± 0.09scaic ± 0.002„, ± 0.04^^,c 

± 0.04„, ± Om^^^Mz) ± 0.002ai ± O.OIA2 1 GeV^ = (3.38 ± 0.11) GeV^ . 



(66) 



The central values are obtained for the matching scale jjo = 120 GeV and the low-energy scale 
fib = 5 GeV. The uncertainty from missing higher order perturbative corrections have been 
estimated by increasing and decreasing the scales /io and fib by factors of 2. Uncertainties 
induced by m^, nib, '^c, C, as{Mz), Ai and A2 are obtained by varying the various inputs 
within the errors given in Table [TJ We assume the errors on C and rric to be fully correlated. 
The total errors are again obtained by adding the individual ones in quadrature. In order to 
show the stability of the zero under change of the b quark mass scheme we collect in Table [2] 
the results we obtained in the IS, MS and pole schemes. 

The total errors are 3.4% and 3.3% respectively and therefore quite small. It is often argued 
that especially the small fi dependence at the zero is an accident and should be increased by 
hand. We argue in the following that the small fi dependence is a reasonable reflection of the 
perturbative error. 

One test of our estimation of the perturbative error on the zero consists in comparing the 
central values and fib dependences of lower-order predictions. The comparison between the 
NNLO+QED, NNLO and NLO extraction of the zero reads (here we consider only the muon 
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NNLO vs NLO 




Figure 2: yUf,-dependence of the forward backward asymmetry for the muonic final state. The 
hnes are the NNLO QCD result; the dashed line corresponds to yU^ = 5 GeV, and the solid 
lines to /ib = 2.5, 10 GeV. The shaded area is the region spanned by the NLO asymmetry for 
2.5 GeV < I2b<l0 GeV. 



channel and quote only the scale uncertainty): 



(3.50 ± 0.10) GeV^ NNLO + QED 
go = <( (3.45 ±0.11) GeV^ NNLO (67) 
(3.11 ± 0.39) GeV^ NLO . 



In Fig. [21, we plot the fib dependence of the FBA in the low-g^ region and compare NNLO 
and NLO QCD results. Here we note that the scale uncertainty is maximal at ~ 1 GeV^, 
decreases smoothly at larger and almost vanishes at the edge of the low-g^ region. Moreover, 
the band between the two solid lines, representing the scale uncertainty of the NNLO FBA, lies 
to a large extent within the shaded area of the NLO scale uncertainty. The same holds true if 
we compared NNLO+QED vs. NLO QCD results. 

In Fig. [H] we show the entire - parametric and perturbative - error band of the full 
NNLO+QED asymmetry over the whole low-s region. The plot shows that the perturbative 
expansion converges nicely everywhere in the low-g^ region. 

The numerical results of Eq. (!67j) and the plots therefore suggest that the variation of the 
scale fib properly describes the size of the missing perturbative corrections. 

As already stressed before, the errors considered here are parametric and perturbative ones 
only; and unknown sub leading nonperturbative corrections of order 0{asA/mb) may give an 
additional uncertainty of order 5%. 



21 



NNLO + QED 
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Figure 3: The full NNLO+QED asymmetry (dashed line) and the total - parametric and 
perturbative - error band (shaded area). 



4.6 New- physics formulae 

Here, we present numerical formulae for the various observables for non-SM values of the high- 
scale Wilson coefficients of the operators P7, Pg, Pq and Piq: 



lohigh 



^high 



2.399 - 0.002576 X(5i?io) - 0.01277 J(5i?7) + 0.0002656 I{5Rj5Rl) 

+0.0004108 X(5i?7(5i?;) +0.002027 X((5i?8) - 0.00003375 J((5i?85^io) 
+0.001676 1{5Rs5Rl) + 0.1079 1{SRg) + 3.022 n{5Rw) + 0.001146 n{5Rio5R;) 
+0.0001236 n{5RioSRl) - 0.0173 n{SRioSRl) - 0.1312 TZiSRr) - 0.0143 n{5Rs) 
+0.0008111 n{SRs5R;) + 0.9375 n{5Rg) - 0.0568 n{5Rg8R*) 
-0.006099 n{5Rg5Rl) + 1.558 \5Rw\'^ + 0.003436 ISRy]'^ + 0.00004162 ISRs]^ 



+0.2231 \5Rg 



X 10" 



(6^ 



2.085 - 0.002576 J(5i?io) - 0.011 J(5i?7) + 0.0002656 liSRjdRl) 
+0.0004108 X{5RrSRl) + 0.002162 I{5Rs) - 0.00003375 I{5RsSRIq) 
+0.001676 X{5R85R;) + 0.09845 I{SRq) + 2.73 n{5Rio) + 0.001146 n{5Rw5R^) 
+0.0001236 n{5Rio5Rl) - 0.0173 ■R{5Riq5RI) - 0.1119 7^(5i^7) - 0.01279 n{5Rs 
+0.0006912 niSRsSP;) + 0.8243 n{5Rg) - 0.04872 n{5Rg5R^) 
-0.005478 n{5Rg5Rl) + 1.411 \5Rw\'^ + 0.002655 \5R7\'^ + 0.00003702 \5Rf,\^ 



+0.2016 \5Rg 



X 10" 



(69) 
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where 



1 + 6R 



7,8 



(oo)ofr 

7,8 



^(00)efr,SM 



-'7,8 



and 



1 + 6R 



9,10 



^9,10 (/^O 
^(ll)SM/ 



(70) 



and we refer the reader to the definition of the Wilson coefficients given in Ref. [30]. 

The analogous formulae for the integrated forward backward asymmetry in bin 1 and bin 2 
are (we give separately the FBA from Eq. fllOj) . A^^, and the normalizing branching ratio from 



Eq. dHD, 

\ '''' /binl 



- 7.877 - 0.84 1{6Rio) + 0.1901 1{6Rs) + 0.1901 1{SRsSRlo) 
+0.02323 J(5i?9) - 8.136 n{5Rw) + 0.02323 X{5R^) n{SRio) 
-6.114 n{5Rio5R*) - 0.6183 n{5RwSRl) + 5.048 n{SRw5Rl) 
-5.942 niSRj) - 0.6059 TZiSRs) + 4.947 n{5Rg) - 0.02323 1{SRio) TZiSRg) 
+0.01569 niSRgSR;) - 0.08693 \5Rio\^ - 0.0128 |5i?9p ] x 10"^ , (71) 



MM/binl 



8.653 + 0.07321 1{6Rj) + 0.0167 1{6Rr6Rl) + 0.002753 1{6Rr6Rl) 
-0.003923 J(5i?8) + 0.01471 X{5Rs5Rl) - 0.03578 X{5R^) + 10.62 n{dRio) 
+0.008387 n{SRio5R;) + 0.001049 n{SRioSRl) - 0.05413 n{5Rw5Rl) 
+0.262 n{5R'j) + 0.02037 7^(5i?8) + 0.05047 n{5Rs5R*) + 2.264 n{5R^) 

10 1 



-0.4778 n{6Rg6R;) - 0.05645 n{6Rg6Rl) + 5.463 |5/2io|^ + 0.2087 |5i?7|2 



+0.002892 l^i^sT + 0.7634 \6Rg 



X 10" 



(72) 



5.459 - 0.7026 J(5i?io) + 0.1426 I{SRs) + 0.1426 liSRgdRlo) 
+0.03706 TiSRg) + 5.472 n{SRw) + 0.03706 X(5i?9) n{5Rw) 
-4.48 n{6Rio6R;) - 0.4645 n{6Rio6R;) + 7.832 n{6Rio6R*g) 
-4.352 7^(5i^7) - 0.4645 7e(5i?8) + 7.622 7^(^i?9) - 0.03706 J((5i?io) 7^((5i?9) 



-0.01172 n{5Rg5R*) - 0.1382 



101 



0.02034 |5i?9p 



X 10 



-8 



(73) 



7.052 + 0.02856 X(5i?7) + 0.005661 1{SR75Rl) + 0.002082 I{SR75R*g) 
+0.01401 J(5i?8) + 0.01112 I{5Rs5R*g) - 0.03512 J(5i?9) + 9.472 7^(5i?lo) 
+0.006338 n{SRwSR*j) - 0.04942 n{SRio5R;) - 0.4596 7^(5i?7) 
-0.05505 7^(5i?8) + 0.01678 TZiSRsSR;) + 2.394 7^(^/^9) 
-0.3527 n{5Rg5Ry) - 0.04016 n{5Rg5Rl) + 4.879 |5i?io|^ 



+0.07016 |5i?7r + 0.6891 \6Rg 



X 10" 



(74) 
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The integrated asymmetries in bin 1, bin 2, and in the whole low-g^ region are then according 
to Eq. (HM: 



( AFB\ _ \ ) binl (nr\ 

V-^/^A^Jbinl - {B....\. . ' 



MM /binl 



(jFB\ _ V^^^^^ )un2 

\ ^'^')\o^ (B ^ +{B ^ ' ^ ' 



5 Summary and Outlook 

In this paper we extend the analysis of log-enhanced QED corrections initiated in Ref. [30] to 
the B — > Xsii decay width in the high-g^ region and to the forward backward asymmetry. 

We give a complete phenomenological analysis of all relevant quantities related to these 
observables, including the ratio TZ{sq) recently proposed in Refs. [24,49]. We propose a new 
approach to the extraction of the zero of the FBA and argue that the scale dependence obtained 
by this procedure is a reasonable reflection of the perturbative error. 

It is well-known that the measurements of those quantities, in addition to the B Xgj 
branching ratio, will allow to fix magnitude and sign of all relevant Wilson coefficients in the 
SM. Since at the end of the current B factories quantities integrated over certain g^-bins will 
be already accessible with high precision, the following observation is important, namely that 
the double differential decay width decomposed according to {z = cos 9) 



dq'^ dz 8 
where 

1^ = Hriq') + H,{q') , ^ = 3/4H^{q') , (79) 

includes a third quantity which depends on a different combination of Wilson coefficients [50]. 
The NNLO+QED analysis of this observable will be published in a forthcoming paper. 
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